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Abstract 

We consider the nonlinear Klein-Gordon equation in R"*. We call multi- 
solitary waves a solution behaving at large time as a sum of boosted 
standing waves. Our main result is the existence of such multi-solitary 
waves, provided the composing boosted standing waves are stable. It is 
obtained by solving the equation backward in time around a sequence of 
approximate multi-solitary waves and showing convergence to a solution 
with the desired property. The main ingredients of the proof are finite 
speed of propagation, variational characterizations of the profiles, modu- 
lation theory and energy estimates. 
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1 Introduction 
1.1 The main result 

We consider the following nonlinear Klein-Gordon equation 

utt - + mu - \u\P~^u = (nlkg) 

where u : M x M'' — > C, to e (0, +oo), and the nonlinearity is i7^-subcritical, i.e. 
l<p<l + -^ifd>3orl<p< +0O if rf = 1,2. 

This equation arises in particular in Quantum Physics where it has been 
proposed as a simple model describing a self interacting scalar field. Mathe- 
matically speaking, the Klein-Gordon equation is one of the model dispersive 
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equations. It is a Hamiltonian equation which is invariant under gauge and 
Lorentz transform and in particular it conserves energy, charge and momen- 
tum. Due to the sign of the nonlinearity, the equation is focusing. At the 
balance between dispersion and focusing, we find "truly" nonlinear solutions: 
the stationary/standing/solitary waves. One can refer to the monograph 
by Nakanishi and Schlag for a complete introduction to the mathematical study 



of equations similar to (nlkg), in particular the study of the dynamics of the 



equation around stationary/standing waves. 



A standing wave with frequency a; S R is a solution of ( nlkg ) of the form 



u{t,x) = e (pi^{x). Such solution has the particularity to exist globally and 
to remain localized at any time. In the physics literature this kind of solutions 
are sometimes referred to as Q-balls. A soliton (or solitary wave) with speed 
V £ M.''', frequency a; € R and initial phase and position G R, aig € K'' is a 



boosted standing wave solution of ( mlkg ) . More precisely a soliton is a solution 
of (nlkg) of the form 



e ^ 



is the Lorentz 



y^uj.vix -Vt- Xo), 

where ipi^^y is a profile depending on uj and v, and 7 
factor. 

We shall consider ground states solitons, i.e boosted standing waves with 
ground states profiles (profiles minimizing a certain action functional, see Sec- 
tion ^ for a more precise definition). The orbital stability properties of such 
solitons have been widely studied (sec [ p| p| , ||, H) || ||, |^ |||). Fur 



ther developments on the dynamics of (nlkg) around solitons can be found in 
the above mentioned monograph 1 38 (see also E^, 39 , Eol pTl ) . Regarding sta- 



bility, it is known that the ground state solitons are stable [EM if the parameters 
are within the following open set 



stab 



{uj,9,v,x) £ R^+^'^i |w| < Vm, \v\ < 1, 



1 



1 UJ^ 

< — 

— rf m 



4 
p-i 



(1) 



Note that Ostab is nonempty only Hp < 1 + |, i-e. the nonlinearity is L^- 
subcritical. 

In this paper we address the question whether is it possible to construct 
a multi-soliton solution for (nlkg), i.e. a solution behaving at large time like 
a sum of solitons. Multi-solitons are long time known to exist for integrable 
equations such as the Korteweg-de Vries equation or the ID cubic nonlinear 
Schrodinger equation (see e.g. ||3^, ^). In the recent years, there has been a 
series of works around the existence and dynamical properties of multi-solitons 
for various dispersive equations, for example the generalized Korteweg-de Vries 
equation jP, p|, 27|, p8| ^ ^ or generic nonlinear Schrodinger equations ^ 
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However, to our knowledge the present paper and 
the recent preprint [|0| are the first works dealing with existence of multi-soliton 
type solutions for nonlinear Klein-Gordon equations (see nevertheless for 
related results on the nonlinear wave equation). 

Our goal is to prove the following existence result for multi-solitons. 



2 



Theorem 1. Assume that 1 < p < 1 + ^. For any TV G N, take a set of 
parameters {ijjj,0j,Vj,Xj)j=i_,,,^N C Ogtab cind let (ipj) be the associated ground 
state profiles :— fuij,vj, and {jj) the Lorentz factors :— (1 — \vj\'^)^^ ■ 
Denote the corresponding solitons by 



TZj (t, x) :- 



'-t+il 



'ipj{x- 



Vjt 



Xj). 



Define 

Vi, := min{\vj 



Vk 



LOi, := niax{|wj|, j 



, = 1, . . 
1,...,7V} 



(minimal relative speed) (2) 
(maximal frequency) (3) 



There exists a — a(d,N) > 0, such that if Vj ^ Vk for any j ^ k, then there 
exist To € M and a solution u to (nlkg) existing on \Tq, +oo) and such that for 
all t G [To, +oo) the following estimate holds 



N 



u{t) -Y^T^jit) + dtu{t) -Y.dtTljit) 



N 



< e 



Remark 1 . During the preparation of this paper we have been aware of the 
work [|o) by Cote and Munoz. Our two results are companions in the following 
sense. In the authors used spectral theory and a topological argument to 
prove the existence of multi-solitons made of unstable solitons. To the contrary, 
we use finite speed of propagation, classical modulation theory and energy es- 
timates to obtain the existence of multi-solitary waves based on stable solitons. 
Merging our results together would give the existence of multi-solitons made 
with any kind (stable or unstable) of solitons. 



1.2 Strategy of the proof 

Our strategy for the proof of Theorem |l| is inspired by the works of Martel, Merle 
and Tsai on multi-solitons of Schrodinger equations |^, ^ (see also |^, ^ where 
similar strategies were enforced). 

We start by introducing the mathematical framework in which we are going 
to work in Section ^. After transforming (nlkg) into its Hamiltonian form (^, 
we list the tools which are going to be useful for our purposes: Cauchy Theory in 
X and x H'^~^, Conservation laws. Finite Propagation Speed, standing 
waves, Lorentz transform and finally definitions of solitons and their profiles. 

Then we go on with the core of the proof of Theorem |l|. We consider a se- 
quence of times T" ^ oo, a set of final data u„ = ^ 7?.j(r") and the associated 



solutions (un) of (nlkg) backward in time. The sequence (u„) provides us with 
a sequence of approximate multi-solitons, and we need to prove its convergence 
to a solution of ( ^lkg ) satisfying to the conclusion of Theorem ^ For this 
purpose, we show that each u„ exists backwards in time up to some time Tg in- 
dependent of n and decay uniformly in n to the sum of solitons (Proposition ^ . 
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Eventually a compactness argument (Lemma ^ permits to show that con- 
verges to a multi-soliton of (nlkg) on [To,oo). Most steps are performed in 
Section ^ apart from uniform estimates whose proof needs more preparation. 

As in the existing literature concerning multi-solitons for (nls) , the proof of 
the uniform estimates relies on several ingredients: coercivity of the Hessian of 
the action around each component of the multi-soliton, modulation theory and 
slow variation of localized conservation laws, energy, charge, momenta. 

We study the profiles of the solitons in Sectio n [i[ We characterize the profiles 
variationally using the conserved quantities of ( nlkg| ) (Proposition ||) . These 
characterizations are not surprising, in view of the profiles characterizations for 
Schrodinger equations. Indeed, in both cases the ground state profiles are at the 
mountain pass level, the least energy level and the Nehari level. Nevert heless, 
to our knowledge, such properties were not derived previously for (ntlkg). Our 
proofs are self-contained and do not rely on the (nls) case. 

After obtaining the variational characterizations, we prove a coercivity prop- 
erty (Lemma ^) for the second variation of the action functional around a soliton 
(linearized action) . To this aim, we study the spectrum of the linearized action 
and prove in particular the non degeneracy of the kernel (i.e. the kernel con- 
tains only the eigenvectors generated by the invariances of the equation, see 
Lemma 1^). It is usually a crucial point in these matters and we obtain it with 
the help of the equivalent property in the Schrodinger case. We underline that 
the coercivity properties are related to the fact that our standing waves are 
stable. 

Coercivity of the linearized action is obtained provided orthogonality con- 
ditions hold. This prompt the question of obtaining orthogonality conditions 
around a sum of solitons, which is resolved in Section ^ via modulation theory. 
The modulation result is twofold. First, it shows that close to a sum of solitons 
one can recover orthogonality conditions (see (|35|)) by adjusting the modulation 
parameters phases, translation and scaling. Second, it gives the dynamical laws 
followed by the parameters (see (^). 

Finally, we define cutoff functions around each soliton and use them to lo- 
calize the action around each soliton. We use these localized actions to build a 
global action adapted to the sum of solitons. Several properties are transported 
from the local actions to the global one. In particular, the global action inherits 
from the coercivity (see Lemma |l^ . Due to errors generated by the cutoff it is 
not a conserved quantity, but we can however prove that it is almost conserved 
(i.e. it varies slowly). We use these properties combined with the modulation 
result to prove the uniform estimates. 

This work is organized as follows. In Section H, we set the mathematical 
work context. Section contains the proof of the main result, assuming uniform 
estimates. In Section ET we establish variational characterizations of the profiles 
and use them to prove a coercivity statement for the hessian of the action 
functional related to a soliton. In Section ^ we explain the modulation theory 
in the neighborhood of a sum of solitons. Finally, we put all pieces together in 
Section ^ to prove the uniform estimates. The Appendix contains the proof of a 
compact injection used in Section ^and interactions estimates used in Section ^ 
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2 Mathematical context 



In this section we introduce rigorously all the necessary mat erial f or our study 
and restate our result in the Haniiltonian formulation for (nlkg), which is a 
more suitable formulation for our needs. But before let us precise some nota- 
tions. We denote by L'?(R'^) the standard Lebesgue space and its norm by ||-||^. 
The space L^(R'^) is viewed as a real Hilbert space endowed with the scalar 
product 

{u, 17)2 —Ke uvdx. 

The Sobolev spaces H''{M.'^) are endowed with their usual norms 
product space i^(R'') x L^(R'') we use the norm 



For the 



l|2 
\ui\\2 



l|2 



with similar convention for H^{R'^) x L'^{R'^), H" 
sometimes us the following notational shortcut: 



ixH^-^iR"^), etc. We shall 



dv 
dxd 



.J 



Finally, unless otherwise specified the components of a vector W E L'^{M.'^) x 
L2(R'^) will be denoted by Wi and W2. 



The Hamiltonian Formalism for the nonlinear Klein-Gordon equation (nlkg) 
is formulated as follows. For W2) € H^{M.'^) x L'^{M.'^) we define the following 
Hamiltonian (which we will call energy in the sequel) 



E 



1|| l|2 



^IIY7 l|2 I "^11 ||2 1 

— Vui W10 

2" ^"2 ^ 2 " ^"2 p+1 



\Ui 



p+i 
P+r 



Define the matrix J 



1 

^-1 

(ui, U2) (u, ut) solves the following equation 



Then u is a solution of (nlkg) if and only if 



dt 



JE' 



(4) 



From now on we shall work only with the Hamiltonian equation (Q). 

Due to this Hamiltonian formulation the energy is (at least formally) con- 
served. In addition, the invariance of (^) under phase shifts and space transla- 
tions generates two other conservations laws, the charge Q and the (vectorial) 
momentum P, defined in the following way: 



Q 



Im 



UiU2dx, 



P 



Re / WuiU2dx. 
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With our restrictions on the growth of the nonlinearity in Theorem |l| (L^- 
subcritical) , it is well-known that the Cauchy problem for (Q) is globally well- 
posedness in the energy space H^{R'^) x L^{R.'^). More precisely, the following 
well-posedness theory holds. 

Cauchy Theory in x L^. Assume 1 < p < 1 + |- For any initial data 
Uq — {u1,U2) E H^(R'^) X L^(M'') there exists a unique maximal solution of (||) 

U EC {R,H\R'^) X L^{R'^))nC^ (R,i2(R'i) X H-\R'^)) . 

Furthermore, we have the following properties. 

Conservation of energy, charge and momentum; for all t E R, we have 

EiU{t)) - E{Uo), QiU{t)) - Q(C/o), F((7(t)) = P{Uo), 

Global estimate; there exist Co > such that //ixls) < Cq\\Uq\\ ^ ■ 

Uniqueness in light cones; If IJ is another solution to on (0,T) for T > 
with C/(0) = Uq in {x : \x ~ xo\ < T} for some xq E R"^, then U = U on the 
backward light cone {{t,x) E[Q,T)xR'^ ■.\x- xo\<T ^t}. 
Continuous dependency upon the initial data; if (Uq ) C H^{R'^) x L^(M^ 
converges to Uq in H^{R'^) x L'^iR'^), then the associated solutions (C/„) of (j^ 
converge in C (/, x L^(R'*)) for any compact time interval I C R to the 

solution U of with initial data C/(0) = Uq. 

For this set of results, we refer to the classical papers by Ginibre and Velo 
or the recent review in the paper by Killip, Stovall and Visan. For our 
purposes, we will also need a more refined result on local well-posedness in the 
slightly larger space H'^iR'^) x H'^~^{R'^) for some s < 1 (see Lindblad and 
Sogge or Nakamura and Ozawa ||36|). 

Cauchy Theory in i?" x H^~^. Let s > be such that either s > d/2 or 
1/2 < s < d/2 and p < 1 + 2^2^- For any initial data Uq = (^5,^2) E 
H'^W^) X H'-'^iR'') there exists a unique maximal solution of M) 

U EC {{-T^,T*),H'{R'^) X H'-\R'')) . 

Furthermore, we have the continuous dependent upon the initial data; if (Uq) C 
H''(R'') X H'-'^{R'^) converges to Uq in H''{R'^) x H'-'^iR'^) , then the associ- 
ated solutions [Un) of (|) converge to U m C (/,7?"(K'') x H'-^iR"^)) for any 
compact time interval I C (—T^^jT*), where U is the solution to with initial 
data U{0) = Uq. 

A useful consequence of the uniqueness in light cones (Cauchy Theory in 
X L^) is the following finite speed of propagation property. 

Finite Propagation Speed. Let U = (ui, U2 ) E i7i(R'^)xi2(Rd) be a solution 
of (^ on (— oo,r*]. There exists Cq, depending only on \\U{T*)\\^x^j^2 such 
that if there exist < e and M > satisfying 

[ |Vwi(T*)|2 + \u,iT*)\^ + \u2iT*)\^dx < e, 

J\x\>M 
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then for any t G [— oo,T*] we have 

|V?/i(0|' + IwiWP + \u2{t)\^dx < Coe, 



I \x\>2M+{T*-t) 

Proof. Let xm be a cutoff function such that 

, , f 1 for |a;| > 2M , „^ „ Cq 

^*^(")^\0 forN<Af I|Vxm|I^<^. 

Define Ut* ,m ■— U{T*)xm and denote by Um the associated solution of (^. By 
assumption, we have WUt" ^mW^i^l^ — ^ ^.nd by the Cauchy Theory in x 
the solution Um exists on M and verifies for alH e M 

\\UMml.^^.<Coe. 

However, by uniqueness on light cones. U m and U coincide on {(t, x) e M x M*^ : 
\x\ > 2M + {T* - <)}, and for any t e (-00, T*) this inipHes 

\Vui{t)\^ + + \u2it)\^dx < WUMimrnxL- < Coe, 

|a;|>2Af+(r*-t) 

which was the desired conclusion. □ 

Lorentz transform. Among the symmetries of (||), we already mentioned the 
phase shift and translation. We consider now the Lorentzian symmetry, defined 
as follows. Take U(t,x) = {ui,U2){t,x) and w € R'' with \v\ smaller than the 
speed of light for (Q). namely \v\ < 1. The Lorentz transform Jii'vU of U is the 
function of (t,x) defined by 



iJ^,U){t,x) : 
where r and y are defined by 



■«i(t,2/) 
l(u2{T,y) -vWyUi{T, y)) 



T — T{t, x) := j{t ^ V ■ x) = —t — 7(x — vt) ■ V, 

7 

y = y{t, x) := X — Xy + ^{x^ — vt) = a; — + (7 — l){x — vt)^. 

Here, the Lorentz parameter 7 is defined by 

1 

7 



and the subscript v denote the orthogonal projection onto the vectorial line 
generated by w, that is 

X ■ v 

It is simple algebra to verify that is Lorentz invariant, in the sense that if 
J7 is a solution of (^, then so is JCyU. Also note that the Lorentz transform is 
invertible with inverse I£-ri- 
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standing waves. Take w <E K. In the Hamiltonian formulation, a standing 
wave with frequency cj is a solution of (^) of the form U(t,x) — e''^*$(^(x). 
Plugging this ansatz for U into (^), it is easy to see that $^ = (^'^•^) must be 
a critical point of _E + uQ, hence a solution to the stationary elliptic system 

The solutions of this system are clearly of the form (^^^ ) , where (pi^ satisfies 
the scalar equation 

~ Aipi^ + {m~uj'^)(p^~\ip^\P~'^ip^=0. (5) 

Solutions to (H) and their properties are well-known (see [|| ^, |l^, ^ and the 
references therein). For every a; G {—^Jm^^pm) there exists a unique, positive, 
and radial function ip^ G C^(R'^) solution of (||). In addition, the function ip^ is 
exponentially decaying at infinity: for any /i < {m — uj'^) there exists C(/i, u) > 
such that 

|</?^(2:)| < C(^,a;)e-v^l^l for all a; G M''. (6) 
Furthermore, any satisfy the scaling property 

p,^{x)^{m-ij^)^ p(^{m~uj^y x), (7) 

where (p is the unique positive radial solution to —Atp + p — \(p\P~^(p = 0. The 
function ip^ is called ground state. In dimension d > 2, there exist infinitely 
many other solutions to (|5]), called excited states. In the sequel, we shall deal 
only with ground states solutions to (jsj). Indeed, our analysis deeply relies on 
properties of the ground states which do not hold for other solutions, in partic- 
ular the stability of the associated standing waves (see Section ^for details). 

Remark 2. It is interesting to notice that, although the presence of the nonlinear 
term permits the existence of states with negative energy, the standing waves 
have always positive energies. Indeed, a straightforward computation assures 
that for a standing wave U{t,x) — e*"*<I>(^(a;) the corresponding energy is given 

by 



The fact that p > 1 guarantees that E{^t^) > 0. 

Remark 3. The scaling property (^ guarantees that the energy of the ground 
states varies continuously with respect to lo. This fact implies that the multi- 



soliton solutions for (nlkg) behave at large time like a sum of solitons that are 



allowed to have different energies. A straightforward computation indeed gives 

\\pX = im~^^)^^^\ml 

\\VipX = {m-u^) ^<p-i) IIV^II^. 
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Now (p is solution of — + <^ — \(pY = ^ such that, by means of Pohozaev 
identity, 

IIY7~I|2 d{V-^) 11-112 

Merging all this information we get the relation between energy and w given by 
where 



p-1 / p(2-d) + 2+d 4-d(p-l) X 

2(p+ 1) ' ^ ' ) 

P+3 2, 2n 



2(P+1) 



which can be rewritten as 

(to(p - 1) + 2uj'^){m - w^)^^^ 



(p+1) 

The monotonicity of g^uj) when uj belongs to Ostab follows easily. 

Solitons. Starting from a standing wave, one generates a new family of solu- 
tions to simply by boosting them using the Lorentz transform. These new 
solutions are the solitary waves (or simply solitons). Precisely, take a frequency 
\uj\ < y/m, the profile ■— (j^^^^ ) (where ip^ is the ground state of (||)), a 
phase e M, and a speed and a position v, xo £ S.'^ with \v\ < 1. The associated 
soliton is 

e'^'+^'^^^^ix-vt-xo). 
where the new profile ^^j.v is given by 

<i>^Jx)^e-^^-^-^( pM + il-l)x,) \ 

\j{iuj(p^{x + (7 ~ l)xy) - vVifi^ix + (7 - ^)Xv))J 

By direct computation, and provided we have noticed that 

it is not difficult to see that ^uj,v is a critical point of 

S := E + -Q + V P. 

7 

With all these preliminaries out of the way, we can go on with the proof of 
Theorem |l]. 



9 



3 Existence of Multi-Solitons 



This section contains the core of the proof of Theorem |l] assuming uniform 
estimates (Proposition |l|) which are proved in Section ^. 

Assume that p < 1 + |. Take S N, {ujj,9j,Vj,Xj)j^i at C Ostab, the 

associated Hamiltonian profiles (as in (||)), and o;^ as in (||), (||), (7^) the 
Lorentz parameters, (Rj) the corresponding soHtons 



Rj {t, x) -.^e '''i $j (x -Vjt- 



and R the sum of the soUtons : 

N 



R{t,x) :=^i?,(t,a;). 



Reformulated using the Hamiltonian expression (1) of ( |nlkg| ), our goal is to 



prove that there exists a — a{d,N) > 0, such that if Vj ^ Vk for any j ^ k, 
then there exist Tq € R and a solution U to (^) existing on [Tq, +00) and such 
that the following estimate holds for all t E [To, +00) 



We are going to define a sequence of approximate multi-solitons and prove its 
convergence to the desired solution of (^). Take an increasing sequence of time 
T" — >■ +00 and for each n let J7„ be the solution to (^ obtained by solving (||) 
backward in time from T" with final data [/„(r") — R{T^). Our proof will rely 
on two main ingredients. First we have uniform estimates for the sequence of 
approximate multi-solitons. 

Proposition 1 (Uniform Estimates). There exist a = a{d, N) > 0, and To G M 

( independent of n) such that for n large enough the solution Un of (^) with 
Un{T^^) = R{T'^) exists on [To,T"] and satisfies for all t G [To,T"] the estimate 



\\Un{t) - < e-"V^-* (9) 

Proposition ^ establishes that the approximate multi-solitons [/„ all satisfy 
the desired estimate on time intervals of the form [To,T"], with Tq independent 
of n. The proof of Proposition |l| is rather involved and we postpone it to 
Section ^ (useful informations for this proof are derived in Sections ^ and ^ . 

The second ingredient of the proof of Theorem |] is an ^ x L^— compactness 
property of the sequence of initial data of the approximate multi-solitons. 

Lemma 2 (Compactness) . Let Tq be given by Proposition ^. For any e > 
there exists such that for any n large enough [/„ verifies 

f |V[/„,i(To)|2 + |{/„,i(To)|2 + \Un,2{To)\^dx < e. 

J\x\>M^ 
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The argument for the proof of Lemma ^ is different from the Schrodinger 
equation case. Indeed, we benefit with the Klein-Gordon equation of the Finite 
Propagation Speed, which is not the case for Schrodinger equations where one 
has to us virial identities (see e.g. Lemma 2]). 

Proof of Lemma |^ The result is a consequence of the Finite Speed of Propaga- 
tion and the uniform estimates of Proposition Indeed, take e > and let T* 
be such that e""V™-'^>*'^* < |. Then it follows from Proposition |l| that for n 
large enough 

ii[/„(n-i?(mi^,,^. <|. (10) 

By exponential decay of the sum of solitons, there exists AI^ such that 

/ |V(i?(ni)P + |(i?(T*))ip + \{R{T*))2\'dx < |. (11) 
Combining (|l0| ) and (^l]), we get 

|V[/„,i(T*)|2 + |CZ„,i(T*)|2 + \U„^2{T*)\^dx < e. 

\x\>AU 

By Finite Speed of Propagation, this implies 

|VC7„4(To)p + \Uns{To)\^ + \UnATo)\''dx < e. 



>\x\>2M, + {T*-To) 

Setting Me = 2M^ + (T* ~ Tq) finishes the proof. □ 

Proof of Theorem With in hand our sequence of approximate multi-solitons 
satisfying the desired estimate, the only thing left to do is to prove that it 
actually converges to a solution of (jj) satisfying the same estimate (^ . 

First of all, we show the convergence of initial data. Since Un satisfies (^ , 
the sequence C/„(To) is bounded in H'^{R'^) x i^^R^). Therefore there exists 
Uo £ ffi(K'') X L^{R'i) such that U„ Uq weakly in H'^{R'^) x L'^(W^). We are 
going to prove that the previous convergence is strong in i?*(IR'^) x H''~^{W^) 
for any < s < 1. Take e > 0. Using Lemma ||, we infer the existence of 

> such that for n large enough 



V(7„.i(To)|^ + \Un.i{To)r + \Un.2iTo)rdx+ 

\x\>M^ 

'' |VC/o,i(ro)|2 + |[/o,i(ro)|2 + \Uo,2iTo)\^dx < J. (12) 

Define Xe ■ '^'^ ~^ [0,1] a cutoff function such that x^{x) = 1 if |a;| < Mg, 
Xe{x) = if > 2A4, and IjVxelL < 1- We have 

\\Un{To)-Uo\\Hs^H^,-. < \\{Un{To)-Uo)Xs\\H^^H'-^+\\{Un{To)-U,){l-Xs)\\H^^H^^^ 
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From the compactnes^ of the injection ff*(ri) ^ H'^~^{n) when fl is bounded 
and (5 > , we infer that, for n large enough and maybe up to a subsequence, 
we have 

\\{Un{To)-Uo)Xe\\H^^H^-^ < |- 

Moreover, by ( p^ 

||(C/„(ro) - Uo){l - Xe)\\H^^H^-. < \\{Un{To) ~ Uo){l - Xe)\\m^L^ < | 

Combining the last three equations gives us 

\\U^iTo)-Uo\\Hs^H^-i <e. 

Hence C/„(ro) converges strongly to in WjM'^) x i/''-i(R'^). 

Let us now show that the solution U of (g) in H^iW^) x i^d^d) ^jth data 
U{To) — Uq satisfies the required estimate. By Local Cauchy Theory of (H) in 
iJ*(R'') X iJ*^^(R''), we have the strong convergence 

Un{t)-^U{t) in i7"(M'') X i7"-i(R''). 

for any t € [To,+oo). In addition, by uniqueness of the limit and since Unit) 
is bounded in 7J^(R'^) x L^(R'') (by (^), we have the weak convergence for 
tG [To, +00) 

Unit) Uit) in H^{«^) X L^iW^). 
Therefore, by weak lower semi-continuity of the x L'^-norm and (0), we have 



\\Uit) - < lim|^fr«W - Rimm^L^ < e-V^''**, 

which concludes the proof of Theorem ^ □ 



4 Properties of the profiles 



Since we will be working mainly within the Hamiltonian formulation of ( nlkg ) , 
it will be convenient to characterize the soliton profiles using the conserved 
quantities. We already mentioned that the profile ^u:,v is a critical point of the 
functional action 



S:=E 

or more explicitly a solution to 



7 



-Q + vP, 



-Awi + mwi — \wi\^ ^Wi + i—W2 

7 

W2 — l — Wi + V ■ VWi 

7 



V ■ V?i;2 = 0, 
0. 



(13) 



proof of this fact is included in the Appendix, Lemma [l7| for the reader's convenience. 
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In this section, we are going to give some variational characterizations of ^cu.v 
and study the Hessian S"{^uj,v)- 

As far as we know, the variational characterizations given in the following 
Proposition]^ were never derived before, although they are expected in view of 
what happens in the scalar setting. The ideas on the relationships between dif- 
ferent variational characterizations used further in this section were introduced 
by Jeanjean and Tanaka in |l^ (see also |Q for related results). We believe 
that these variational characterizations of the profile are of independent 
interest. 

For the purpose of constructing multi-solitons, the main result of this section 
is the coercivity property given in Lemma ^j. The proof of this result relies on the 
variational characterization of the profiles as well as on their non-degeneracy, 
which is given by Lemma |7[ We shall follow closely the presentation made in 
for the standing waves of nls. 

4.1 Variational Characterizations 

We define the mountain pass level by 

MP ;= inf sup S{r]{s)) 

where F is the set of admissible paths 

F := {ry e C ((0, l),i/i(M'^) x L^R^)) ; 77(0) = 0, S{v{l)) < 0}. 
We define the Nehari constraint for W G H^iR'^) x L'^{R'^) by 

I{W) := {S'{W),W) 

and the Nehari level by 

NL mm{S{W); I{W) ^0,W ^ 0}. 

We also define the least energy level by 

LE := min{5(M^); W e H^R'^) x L^iR"^), W ^ 0, S'{W) = 0}. 

Proposition 3. The profile ^i^^y admits the following variational characteriza- 
tions: 

^(^^ „) = MP = NL = LE. 

Let us start by proving using mountain pass arguments that S admits a 
critical point. Then we will show that this critical point is at the mountain 
pass level and also at the least energy level and at the Nehari level and we will 
identify it with $^^^1,. 

Lemma 4. There exists ^ G H^iR'^) x L'^{R'^) a non-trivial critical point of S, 
i.e. 

*7^0, S"(*)=0. 
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Before going further, we make the following useful observations on the for- 
mulation of S: for W = (wi, W2) G H^{R'^) x L^(R^) it is simple algebra to see 
that 

SiW) = h\Vwi\\l-hvVwi\\l + \ (m-'^) \\wi\\l~^-vlm f wiVw^ 
2 2 2 \ Y J 2jJ^d 

_L II Y7 I l|2 1 II llP+l 

2 " 7 "-^ p + 1 

We can remark further that if wi is such that wi{x) = e~'^'^'^^'^wi{x+ {-f— 
then we have 

SiW) = ^(^^\\Vw^\\l + ^{m,-u;^)\\w^\\i^ 



-\\v ■ Vwi - i-wi + W2\\2 - ^;^ll^"illp+r (14) 



We shall also use the following Lemma at several occasions. 

Lemma 5. For any e > there exists S > such that for any W — (wi, W2) G 
H^IR'') X L2(Rrf) we have 

e\\wi\\l. + \\v ■ Vw^ - I'^w, + wAl > S\\W\\l,^^.. 
Proof. We only have to make ||t«2|l2 appear. Write 



W2 = a [ V ■ Vwi — i—wi j + W2 , 



where a e M and ■ Vwi — i^wi, = 0. We have 

II?; • Vwi — i—wi + W2II2 = (1 + Q^)^||^^ ■ Vwi — 'fi'iHj + Il''^2"ll2 (1^) 

II l|2 211 V7 ||2 , II _L||2 

IIW2II2 = " 11^^ • VWi - i — + \\W2 \\2- 

There is a possible degeneracy in (^5|) if a = — 1. but we can compensate it by 
using a piece of 

^ II n2 II „ ||2 ^ 

2 11^1 11^1 + \\v- VWi - i—wi + W2II2 > 

C ((1 + af + I) ||« • Vwi - i^will' + C\\w^\\l > C\\w2\\l. (16) 
The desired inequality is then a direct consequence of (Kw. □ 
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Proof of Lemma ^j. Step 1 : Mountain-Pass geometry. 

We claim that the functional S has a mountain-pass geometry, i.e 

MP > 0. 

We start by showing that r is not empty. Indeed, take W = (^1,^2) G H^{M."-)x 
L2(R'^) and s > 0. Then using ^ we see that 

"2 ^1/1„_. „, 1 



2 , (' 2^ II ~ ii2 



111 

+ — \\v ■ Vwi — i—wi + W2\\r, 7 ^"1 



2ir ' " "'-^ " "^^112 / ;r+I"""'p+i 

Therefore if s is large enough we have S{sW) < 0, hence the path s H> S{-^W) 
belongs to F provided C has been chosen large enough. 

To show that MP > 0, it is enough to prove that there exists a function 
/ : R+ ^ M such that f{s) > for s close to and S{W) > VFH^i ^^2). 
Using ([1^), the continuity of wi — > wi in H^{]&.'^) and Sobolev embeddings, it 
is easy to see that there exists e > such that 

S{W) > UwiWl, + ■ Vti'i - i-wi + W2\\l - C\\wi\C^\ 
Z Z 7 

From Lemma [s] we infer that there exists S > such that 

siw)>~s\\w\\l,,^.-c\\wC,',^,. 

This imphes that S{W) > if ||W^||/jixl2 small and there exist C > 0, 5 > 
such that S{W) > C > for = ^- This implies MP > and has 

a mountain pass geometry. 
Step 2 : Existence of a Palais-Smale sequence. 

From Ekeland variational principle (see e.g. [^) and Step 1, we infer the 
existence of a Palais-Smale sequence Wn = (wi_„, W2.n) at the level MP, i.e. 

S{Wn) -> MP, S'{Wn) 0, as n +00. (17) 

Step 3 : Non-vanishing of the Palais-Smale sequence. 

Assume by contradiction that the sequence Wn is vanishing, more precisely 
for any i? > we have 

lim sup / {\wi,n\'^ + \w2.n\^) dx ^ Q. 

"^^^"^ ym-^ J\x-y\<R 

Take e > and i? > and let n be large enough so that 



sup / (|u;i,„p + |u;2,nP) rfa; < £. 

yGR-^ J\x-y\<R 
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Recall Lions' Lemma (see [^): for any w G H^{R''-) we have 

\\w\\lX\<c{snv[ \w\^d^ \\w\\\,. (18) 



l\x-y\<R 

Therefore, for n large enough, and using (Q) and Lemma ^, we get 

S{Wr.)>C\\Wn\\l.^L.-e\\w^,r.\\\., 

which implies (if e has been chosen small enough) that [Wn) is bounded in 
H'^iW^) X L'^iW^). This boundedness has two consequences: as n — +00, we 
have 

(^'(PF„),W„)^0 and ^0. (19) 

where the first limit is due to the fact that [Wn) is a Palais-Smale sequence 
(see (^7|)) and the second limit comes from (|8|). However, we have 

and therefore (^91) implies 



lim S{Wn) = 0, 

n— ^+oo 

which enters in contradiction with lim„^+oo S{Wn) — MP > 0. Therefore the 
sequence {Wn) is non- vanishing. 

Step 4 ■ Convergence to a critical point 

Since Wn is non-vanishing, there exists R,S > and (y„) C M'* such that 
for n large enough 

/ \wi,n\''dx>S. (20) 

If we substitute Wn{- — Un) to Wn (keeping the same notation), the sequence 
(Wn) is still a Palais-Smale sequence and keeps the same properties. In particu- 
lar, as known from Step 3, (Wn) is bounded in H^{M.'^) x L'^{R'^), and therefore 
we have the existence of * G i?i(M'') x ^^(R'*) such that W„ ^ * weakly in 
i?^(R'') X L'^{W^). Since Wn is a Palais-Smale sequence and S' is continuous, we 
have 5"(5') = 0. Hence we only have to show that ^ is non-trivial. This is a di- 
rect consequence of (20) and the compact injection -ff^(|a;| < R) ^ L'^(\^ < ^)- 
Hence is a non-trivial critical point of S and the proof of Lemma ^ is fin- 
ished. □ 

We turn now to the variational characterizations of the critical point ob- 
tained in Lemma ^. 

Lemma 6. Take 'J the critical point of S found in Lemma Q. The following 
equality is satisfied. 

5(*) = MP = NL^ LE. 
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Proof. Let us start by showing 

S'(*) < MP. (21) 

Using S"(v[') = 0, wc have 
p+l 

+ \\v ■\/wi-i'^wi+W2\\fj < ■■. (22) 

Recall that wi is such that wi{x) = e~'"'''"'^?2)i(x+ (7— l)xy) (see (p^). Using 
the weak lower semi-continuity of the norm, we can continue the inequaHty 
started in (^) by 

II ^ -"^ l|2 

+ \\V ■ VWi^n - l — Wi^n + W2,n\\^ 
7 

= hminf (s{Wr,) - (5'(W^„), W^„)) • (23) 

Since {Wn) is a Palais-Smale sequence we have 

lim (s{Wn) - {S\Wn), Wn)] = MP, 

and we can conclude from (|2|) and (|3|) that * verifies (|l]). 
We continue by showing that 

MP < NL. (24) 

Take an element of the Nehari manifold W E H^iR'^) x L'^{M.'^), I{W) = 0. 
The idea, as in [|l8[ |l^, is to construct a path in F so that S{ri{s)) achieves its 
maximum when ri{s) — W. It is easy to see that for C large enough the path 
rjc defined by t]c{s) = CsW fulfills our needs. Indeed, we have 



d_ 

ds 



S{sW) = s(-(||Vwi||2 + ("I - w') ll^ill 



I I .UJ I I 2 Y) 111 I I M + 1 

+ • Vwi - «— wi + W2II2 - s ll"'illp+i 

In particular, §-^S{sW)\s=i = I{W) = 0. Therefore ^S'(sll^) > for s e (0, 1) 
and -^S{sW) < for s > 1. Hence the path S{r]c{s)) achieves its maximum 
when s — ^ and 77(^) = W. Therefore, 

ML < S{W), 
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and since this is true for any W on the Nehari manifold this proves (p4|). 
It is easy to see that 

NL < LE. (25) 



Indeed, any solution W of (i.e. any critical point of S) satisfies the Nehari 
identity I{W) = 0. Thus the infimum for NL is taken on a larger set than the 
infiniuni for LE, hence (|25|). 

Finally, as a direct consequence of 5'(5') ~ and the definition of LE we 
have 

LE<S{^). (26) 
Combining (|2l|), (H), (||), finishes the proof of Lemma | □ 

Proof of Proposition |^. In view of Lemmas ^ and ^, the only thing left to prove 
is that ^ = ^u:,v Let us first see the case v = 0. Since 5* = (-01, ^^2) is a critical 
point of S, we have ?/'2 = iujij^i. Therefore, since ip^ is a ground state of (|^), we 
have 

p+i 



5(vl') = 2 llV^ill^ + - (m - ll^ill^ - U41X\ 



> 



Therefore when w = 0, we indeed have = $cj.o- Let us now treat the case 
V 7^ 0. Let ifji be such that ipi{x) = e~*"^"'^?/'i(a; + (7 — l)a;.„) and define 
■ip2 ■— iwV'i. Then ^ := (V'i,V'2) is a solution to ( p^ ) with v = 0. Indeed, it is 
not hard to see that 

-AV^i + i^V2 - • VV'2 = e^*"''^-^ (-AV^i - w^i/ii) . 

Hence 

^ + ,^g)(^r) > i^E + c^g)($^.o) = 5($^,„). 

7 7 

This implies that = <I>^.„ for any v and finishes the proof of Proposition ^. □ 
4.2 Kernel 

Lemma 7. T/ie following description holds for the kernel of S"{^aj,v)- 

Ker(S"' ($„,„)) = Span{i$„^„, V$„,„}. 

Proof. The inclusion D is easy to obtain. Indeed, due to invariance by transla- 
tion and phase shifts, for any G R and y ^ M.'^ we have 

The result is obtained by deriving with respect to 6 and y at 6* = 0, y = 0. The 
reverse inclusion is much more delicate. We shall rely on existing results for 
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standing waves of NLS to prove it. First remark that if = (ifi, W2) belongs 
to the kernel of S"{^uj.v), then it satisfies 

+ \^lv\^~^Wi + i-W2 - V ■ VW2 ^ 0, 

7 

W2 — i—Wi + V ■ Vwi = 0. 

7 

Here, we have denoted by ,„ the first component of ^u,v, i-e. 
Take W := {wi, W2) such that 



W1{X + (7 - l)Xy) 

jW2{x + (7 - l)xy) - ^vVwi{x + (7 - l)Xy) 



It is a lengthy but straightforward computation to verify that W satisfies 

W2 — i(^wi = 0. 

Remembering now that — ip^^ and using the second equation to substitute 
in the first we get 



W2 — lUJWi. 



(27) 



Fortunately we arrive on a known ground: it is well-known since the celebrated 
work of Weinstein and Kwong (see |^ for a modern short proof of this 
result) that the only solutions to (^7[) are 



W2 



€ Span 



Coming back into the original variables, this implies that 

W e Span{V$^,t,;i$^,i,} 

and finishes the proof. 



Span{V$^,o;i$;^,o}- 



□ 



4.3 Coercivity 

The proof of our result relies on the fact that the solitary waves we are consid- 
ering are stable. In particular, we have at our disposal a coercivity property on 
the Hessian of the action S related to the soliton profile ^^u.v which allows us 
to control the difference between a soliton and a function in a neighborhood of 
its orbit. The coercivity property is the following. 
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Lemma 8 (Coercivity). Assume p < 1 + j and let uj <E ^ and w G K'* be such 

that j < ^ < 1, and \v\ < 1 (i.e. they are compatible with Ostab defined 

in (jlj) ) and let be the associated ground state. There exists S such that for 
any W € H^(W'') x L^(K'^) satisfying the orthogonality conditions 

(W, V$^,,)2 = {W, I J$..,.)2 = {W, = (28) 

we have 

H^Aw)>s\\w\\l^,^,. 

where for brevity in notation we defined 

H{W) := (5"($„,,)VF,1¥). 

Similar results date back to the work of Weinstein |T) for nls equations. 
These ideas were later generalized by Grillakis. Shatah and Strauss in a 
abstract setting. More recently, Stuart |Q described precisely the orbital sta- 
bility of solitons of NLKG using also a coercivity statement, but with different 
orthogonality conditions and a slightly more complicated proof. 

Proof of Lemma^ Step 1 : Analysis of the spectrum of S"{^uj,v)- 

We first remark that, due the exponential localization of ^ui,v, the operator 
S"{^uj.v) is a compact perturbation of the self-adjoint operator 

By Weyl's Theorem, S"{^u:,v) and L share the same essential spectrum, that 
we now analyze. Observe that for W = (wi,W2) S H'^{M.'^) x H^{M.'^) we have 

iLW.W) ^ llVunllUnMl + lW.f, - 2^ + 2„ ■ (V„„„„, , 

which, after some factorizations (similar to those used in (p^)), we can rewrite 
{LW,W) = - (||Vi&i||2 + (to- w^) llwills) + \\v ■ Vwi - i-wi+W2\\l 



where wi is such that wi{x) — e *"^"'^wi(x + (7 — l)xv). From Lemma ^, we 
see that there exists S > such that for any W € H^{R'^) x H^{M.'^) we have 

{LW,W)>5\\W\\l,,^.. 

This implies that the essential spectrum of S"{^i^A is positive and away from 
0. The rest of its spectrum consists in a finite number of isolated eigenvalues. It 
turns out that from the variational characterization of ^u,vi we can infer that 
S"{^uj,v) has Morse Index 1, i.e. it admits only one negative simple eigenvalue 
(see e.g. [Q). We denote this eigenvalue by —A < 0, and ^! an associated 
normalized eigenvector, i.e. S"'($tj,ti)^ = — A^t and ||^||i2xL2 = 1- 
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Step 2 : A positivity property. We prove now that ifW€ H^{'R'^) x L^(R'^) 
satisfies the orthogonality conditions (pq), then 



{S"{^^,,)W,W) >0. 

A particular vector associated with S"'($tj,„) is A(^$(^_„, where A^j := 
Indeed, deriving S'{^i^^y) = with respect to uj we get 

5"($..,,)A„$^.„ = --Q'($...„) = --zJ$„,„. (29) 

7 ■ 7 

This impHes, using 

7 ' ■ 7 ■ 



UJ ,|2\ (to — 1^ 



2)^-f / 2^2 + I 



= A^ i^^iMt) = "'^ — I — r " - 1 1 ii^ii' < (30) 

where the last inequality follows from the fact that to is compatible with Ostab 
(see the definition (|l|)). It is easy to verify that A^^^^^^y is orthogonal to the 
kernel of S"{^uj.v), namely that we have 

(A<^$^,i,,i$„,^,)2 = (A^$^,.„, V$^,„)2 = 0. 
Let us write the orthogonal decomposition of A^^i^ ^ along the spectrum of 

A^<^>^^y = a* + n, (31) 
where a ^ and 11 is in the positive eigenspace of S"{^i^^y), in particular 

(y'(<i>^,„)n,n) >5||n||^,,^.. 

From (|0l) and (|l|), we infer that 

- Xa^ + (5"($^,„)n, n) =. (5"($„,„)Ac.$c..., A^$^,„) < (32) 

Take now W G H'^{R'^) x L^iR'^) satisfying the orthogonahty conditions (||). 
We also write the orthogonal decomposition of W along the spectrum of S"'($tj^^): 

W = IS^i + E, (33) 

where /3 € R and S is in the positive eigenspace of S"(^ui.v)- If /3 = 0, the 
conclusion follows, so we assume (3^0. Using (^8|), (|29|), ( pll) and (|33|), we 
have 

= (5"($^,„)A^$^,., M^) = -Aa/3 + ($„,„)n, S) 

Note that on the positive spectral subspace of 5"($^_i,), Cauchy-Schwartz in- 
equality holds: 

(5"($„,„)n,s)' < (5"($„,„)n,n) (5"($^,„)5,5) . 
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Therefore, 

{S"{^^,,)W,W) = -A/32 ^ > 

+ (5"($„,.)n,n) +^^-° 

Step 3. The coercivity property. 

Assume by contradiction that there exists {Wn = {w",W2)) C H^{M.'^) x 
L'^(W^) satisfying the orthogonahty conditions ( p^ ) and such that 

\\Wn\\m^L^='^ and hm (5"($^,„)VK„, W„) = 0. 



Recall that, as for (14), we have 



(5"($„^,)W^„, VF„) = i (llV^II^ + (m - C^2) ||^n||2^^||^.^^^_^^^^^^^||2 

where is such that w"(a;) = e~"^''''''^w"(a; + (7 — l)a;«)- 

Since is bounded in ^^(M'') x L'^(W^), there exists e i/i(M'') x 

such that 

Wn^W a.sn^ +00 weakly in H^{W^) x i^^Rd), 
On one hand W must satisfy (^) and from Step 2 we have, iiW^O. 

{S"{^u,,,)W,W) > 0. 
On the other hand, by weak convergence and exponential decay of ^uj,v we have 
{S"{'i>^.^)W, W) < liminf (5"($^,,)W„, W^) = 0. 

n— ^+00 

Therefore must he W = 0. However, in this case it would implies 

and since ||W^n||^fixL2 = 1, we would have (using Lemma |[) 

(^"($^,„)W„,VK„)>5>0, 
which is a contradiction. □ 
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5 Modulation Theory 

The use of a coercivity property similar to Lemma ^ but adapted to a multi- 
soliton (see Lemma |l^) will require to deal with orthogonality conditions. These 
orthogonality conditions will obtained by modulation. 

Given parameters (e, L), consider a neighborhood of the sum of solitons 



U{e,L) ■={U e H^{R'')xL 



inf 

3 = 1,. ..,N 



N 



< e 



The main result of this section is the following. 

Proposition 9 (Dynamical Modulation). There exists e, L,C,C > such that 

for any < e < s and L > L the following property is verified. 

Let U{t,x) — {ui,U2){t,x) be a solution of (^) satisfying on a time interval I 

UeU{e,L), for all tel. 

For j = 1, . . . , N , there exist (unique) functions 

Bj-.I-^R, UJ-j : I {-y/m,^/m), Xj-.I-^R'^, 

such that if we define Rj(t) and T{t) by 

N 

i?,(t) = e'^^^W$^^.(t),„^.(.-i,(i)), T{t)^U{t)-J2R,{t), (34) 



then T satisfies for all t e I the orthogonality conditions 



Moreover, for all t e I we have 

N 



-ujj\< Ce, 



0, i = 1, . 



L 



,iV. 



(35) 



Xj+i - Xj > -, j = l,...,N -I, 



and the derivatives in time verify 



N I 



l3 



(36) 



The proof of Proposition || relies on the following Lemma. Note that this 
lemma is valid for time-independent functions. 
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Lemma 10 (Static Modulation). There exist L,C,e > such that for any 
L > L, < e < e, the following property is verified. 
For j = 1, . . . , N , there exist (unique) functions 



:U{e,L) 



ujj\U(£,L)^{—^/m,^/m), Xj:U{e,L) 



such that if we define and T by 



N 



then T satisfies the orthogonality conditions 
Moreover, 



0, j = l,...,N (37) 



N 



I'^IIhixz,^ +y]l'^j -'^jI ^ C-e, Xj+i-Xj > ^, j = 1,...,N -1. 

i=l 



In the proofs, we will use the notation A<^^. for the scaling operator, i.e. 

Proof. We start by proving the lemma in a ball. Take £ > 0, L > 0, (t^j)j=i,...,jv C 
M and (Ci)i=i,...,jv C M'^ such that ^^+1 > + L hv j = 1, . . . ,N - 1. Let B(£) 
denote the ball of H'^{R'^) x L'^{R'^) defined by 



B(£) = lUeH^iR'') X 



U-Y^e'^^^jix-^,) 



< £ 



Define po := (i?i, . • . , ^Jw, wi, . . . , w^, 6, • • • , ^iv), and let «P C x x (M<*)^, 
be a neighborhood of po- We denote by 

p = (6'i,...,6'jv,n7i,...,rojv,yi,...,2/jv) 
a generic element of «p. We define the functional F:^x B{e) (M(''+2))^ by 

F{p,U) := I FkMP.U) 
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where for fc = 1 , . . . , we have set 

N 



2 



Here, is the translation by y, i.e. Tyv(x) — v{x — y). We clearly have 

The lemma inside the ball will follow from the Implicit Function Theorem if we 
prove that 

' p = po , [/ = ^ e'^^ $j {x - ) is invertible. (38) 



dp 



The computation of the derivative is not very hard. Many terms will be made 
small using the exponential decay of the profiles. Other will cancel due to 
orthogonality. We will essentially be left with a diagonal matrix with nonzero 
entries, hence the invertibility. We give only some representative calculations. 
Let's start by 



09, 



3 = 1 



When j — k, we readily have 



AT 



Il2- 



i=i 

Assume now j ^ k. Then, by exponential decay (see (||)), we have 
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Therefore, since \^j — £^k\ > L, this impHes 



J ^ j=l 



< Ce-- J ^. (39) 



This quantity can be made as small as we need by increasing the value of L. 
For the derivative with respect to wj, we have 



When j ^ k, this quantity can be made small as in (^9|). For j = k, since 
ipuj G M, we simply have 



dwk 



\ 1 

.7 = 1 ^ ~ 



All other computations follow from similar arguments and we finally find that 



where DIAG is a diagonal matrix with nonzero entries on the diagonal. There- 
fore, for L large enough, we have the desired invertibility property ( |3^) and 
the Implicit Function Theorem implies the result inside the ball. Since any 
U ^U{e, L) belongs to some ball B{e), the existence part follows in the cylinder 
U{e,L). To show uniqueness, one has to prove that the functions obtained are 
independent of the ball chosen, we leave the details of this argument to the 
reader. □ 

Proof of Proposition |^. The first part of the statement follows from Lemma LO] 
(except the regularity that follows from other regularization arguments, see |p9|), 
hence the main thing to check is (^6|). We first write the equation verified by 
T. Recall that U satisfies dtU = JE'{U). We replace U by X)^! Rjit) + 
in the previous equation to get 

N , N 

such that it follows 



£T + m{T) + o(e-3"\/™-'^^**) (40) 
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where £ is the Hnearized operator defined by 

A + m 0\ f{p- l)\R^\P-^Rj Re{Rf) + \Rj\P-^ 





and *Tt(T) is the remaining nonhnear part. To write this equation, we have used 
Lemma the fact that 



N \ N 



E' [^Rj] ^^E'{Rj) + o(e-3"%A^^*t)^ 



AT 



and that Rj is a critical point of _B + + - P. An analogous computation 
is derived in all details in Lemma 

Take now the scalar product of (^o|) with iJRk- By using Lemma |8|, the 
definition of Rk , and the orthogonality conditions (^) it follows that 

(vRk,iJRk)^ - (tRk,iJRk)^=0 

(Kc^^Rj.iJRk)^ = {vkj.iJRk)^ = (iR,,iJRk)^ = 0(g-3ay^;r^..t) ^. ^ ^ 
(5tT,zJ7?fe)_^ = -(T,9tzJi?fc)^. 
Therefore 

+ 0(||T||^,,^.)+0(e-3"V^-*), (41) 
where the term hu^QiRk) comes from 

(^Ac^Rk,iJRk)^ = (^Rk,iJAc^Rk)^ ^ Acj^Q{Rk). 
Note that by exponential localization 

(£T,iJi?fe)_^ = (T,£*a7?fe)_^ = (t, {JE"{Rk))*iJR}}j ^ + 0(e-3"V"--'^**) 

We want to use the fact that belongs to the kernel of S"{Rk) (see Lemma 
namely that 

J (^P" + ^Q" + ffe • P"^ (Rk)^ iJRk = 0. 
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To this aim, we use the definition ( |34D of Rk, to compute the following time 
derivative and make the missing parts appear. 



T, dtiJRk] = — fx, {JQ"{Rk)riJRk) + vk ■ fx, {JP"{Rk)riJRk 

^ - dA 
Ik 



T,JRk) +ivk-dtik)-{T,iJVRk 



(T,iJAc,,Rk)_ 



Therefore, (ETl) gives 



dtCjk[AoMRk)) - {T,JRk)^ + {vk-dtik)[T,iJVR, 

+ dtCok (T,tJAo,Rk) ^ + 0{\\T\\l,^^,) + 0{e-'^V^-*t). (42) 

Take the scalar product of (^) with -^Rk for = 1, . . . , d to get from similar 
arguments 



{vk - dtik)^ 



+ 0(||T||^,,^.) + 0(e-3"V^-*). (43) 



Conversely to what happened for (^), we do not expect to have a cancellation 
on the linear term. We just estimate it by 



Therefore ([i3[) gives 



{vk - dtXk)- 
a 



^Rk 



= 0(||T||^,,^.)+0(e-3"V^''**). 



(44) 



Finally, take the scalar product of ( |40| ) with iRk and argue as previously to 
obtain 

(5A-^fe)||i?fc|l2 = 0(||T||^,,^.) + 0(e-3"V^-**). (45) 



Putting together (|42|), (^4[) and (^5|) we obtain a differential system for the 
modulation equations vector Mod{t) :— {dtUJj, dtOj — ujj^dtXj — Vj)j=i,...,N of 
the form 

A ■ Mod{t) = M(T) + +0(e-3"V™-"^**), 

where |Af(T)| < C||T||^i^^2- As long as the modulation parameter do not 
vary too much and ||T||^i remains small, A is invertible (it is of the form 
DI AG + small with DIAG a diagonal nondegenerate matrix) and we can deduce 
that 

\Mod{t)\ < C||T||^i^^. + o(e-3"V^^-**). (46) 
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Coming back now to (Q), it is now easy to see that in fact we can improve in 
part the previous estimate into 

N 

J2\d,Co,\ < C||T||^,,^. +0(e-3"v/^-*). (47) 

This improvement is due to our choice of orthogonahty conditions. Combin- 
ing (^) and ( p7| ) gives the desired result. □ 



6 Uniform Estimates 

This Section is devoted to the proof of Proposition |l|. We essentially follow the 
same line as in the Schrodinger case pO| . Since our approximate multi-solitons 
have final data Un{Tn) — R{Tn), they satisfy the desired estimate at least on 
some interval [T„ — (5, T„]. Thus the idea is to reduce things to a bootstrap 
argument: Proposition |l] is a consequence of the following proposition. 

Proposition 11 (Bootstrap). There exist a ~ a{d,N) > 0, and Tq G M (in- 
dependent of n) such that for n large enough the following bootstrap property 
holds. For t^ G [Tq, T"], if [/„ satisfies for all t £ [t\ T"] the estimate 



|!C/„(t) - i?(t)||^.,^. < e-"V^'^**, (48) 
then it will also satisfies for all t £ [t^ , T"] the better estimate 



||C/„(t) - Rmm.L^ < ie-V^''*'. (49) 

Let us quickly indicate how to obtain Proposition ^ from Proposition [Tl| . 

Proof of Proposition |^. Proposition pi] i mplies Proposition ^by means of a clas- 
sical continuity argument (see e.g ||30|). First, let us notice that the map 
1 1-^ Unit) £ H'^{R'^) X L'^iR'^) is continuous. Second, let us define 

t* := inf{T e [To,r"] such that (||) holds for all t £ [r,r"]}. 

Recall that C/„(T") ^ R{T''), therefore we have Tq < t* < T". Our purpose is 



to show that t* — Tq. Let us suppose that t* > Tq. Thanks to (49) we get 



\\U„it)^Rit)\\^,^^,<^e-'^V^-*\ 
for all t £ [t* , T"]. By continuity it exists > such that 

for all t £ [t* — Si,T'"-]. This contradicts the definition of t* and finishes the 
proof. □ 
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Hence now we only have to prove Proposition |Tl|. For the rest of the paper, 
we make the following assumption. 

Bootstrap Assumption. Let Tq > to be determined later and assume that 
there exists t^ G [Tq, T"] such that Un satisfies for all t £ [t\ T"] the estimate 

||C/„(t) - Rmm^L-^ < e-"V^-**. (50) 

We want to prove that in fact (^0|) holds with the better constant ^ on the 
left hand side. 

To prove Proposition |ll|, we need a way to control the difference between 
the sum of solitons and the approximate multi-soliton Un- If there is only one 
soliton. it is known since the ground work of Weinstein [ pO| that the coercivity 
property of the hessian of the action functional (Lemma provides a mean to 
control the difference between a soliton and a solution close to the orbit of the 
soliton. As in ^, ^ ^ , we are going to generalize such a property to the case 
of N solitons. To that purpose, we define localized versions of the conservation 
laws around each solitons and prove that a coercivity property also holds for 
the functional action related to the multi-solitons. 

6.1 Bootstrap 

In this section, we prove Proposition |ri|, assuming three intermediate Lemmas 
proved in the later sections. 

First of all, we begin by selecting a particular direction of propagation. De- 
fine the application J7 : E>'^~-^ — >• R by 

Let ei be such that 

fi(ei) = max {17(e), e e S''"^} > (51) 

Here, the sup is a max since we are maximizing a continuous function on a 
compact set. Let us prove the last inequality. We have 

Each set composing the union on the right is of Lesbegue measure, therefore so 
is fi-i({0}). Hence E,"^-^ \ n-^{{0}) ^ and this proves (|l|). We can complete 
ei into an orthonormal basis (ei , . . . , e^) of M'' and we infer from (^Tj) that there 
exists a > such that for any j ^ k, 

\ivj-Vk,ei)iiid\>a\vj-Vk\- (52) 
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Since is rotation-invariant, we can assume that (ei, . . . , e^) is the canonical 
basis of R'^. Calling vj the first component of the j-th velocity vector, up to 
reindexing the solitons, we can assume that 



v} < vl < ■■ ■ < 



N- 



The localization works as follows. We first define a partition of unity {4>j)j=i^...^N'- 
take ip a cutoff function such that 

ip{s) = for s < -1 and ^(s) = 1 for s > 1, < < 1 in [-1, 1], 



3C > 0, V.S e M, < CV?(s). (53) 

Define 

'/'j = ~ "^j+i for j = 1, . . . , iV - 1 and (/)jv = V'A'- 
We consider the following localized action functional for W = (wi , W2) G ^ (K'*) x 

L2(]Rrf) 

Af JV 

where for j — l,...,N we have defined the localized energies, charges and 
momenta by 

Ej{t,W):^^[ |Vwip0jdx + ^ /" \wi\^(f>jdx + l [ \w2\^(l>jdx 

\wi\P+'^(l)jdx, 



p+l 

Qj{t,W) := Im / wiil!2<pjdx, 
Pj{t,W) ■.= Re VwiW2(l)jdx. 



Since [/„ verifies (50), we can assume that Tq is large enough, so that J7„ 
satisfies the hypotheses of Proposition ^ and thus there exists a modulated sum 
of solitons R = X^jLi ^^^^ T„ verifying the orthogonality conditions ( ^5| ) 
such that 

C/„(t) = i?(i) + T„(i), 



||T„||^,,^. <Ce-"V^''**. (54) 
Let us define the localized linearized action for R by 

3 = i 
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It turns out that T-Ln is inheriting the coercivity property of the hessian of 
the action around a single soliton (Lemma ^) . 

Lemma 12 (Coercivity). There exists C > such that for all t E [To,r"] the 
localized Hessian verifies 

nniTr.it), T„{t))>C\\TrM\m^L^- 

In addition, since Sit, Unit)) is made of localized versions of conserved quan- 
tities, it varies slowly. 

Lemma 13 (Almost conservation). Fort e [t*,T^], we have 
d 



Q-^Sit,Unit)) 



We also have the following Taylor-like expansion for Sit, Unit)). 

Lemma 14 (Taylor-like expansion). The action Sit, Unit)) satisfies for t € 
[t*,T"] 

N , . 

Sit, Unit)) = [E^RM + -QiRji^)) + ^1 ■ PiRjit)) ) 

i=i ^ / 

+ UniTnit), Tnit)) + o(e-2" V^^"** ) , 



With Lemmas [L2, [13 and 14 in hand, we can now conclude the proof of 



Proposition 11 



Proof of Proposition |7J. The first step is to show that 
Indeed, thanks to Lemma [13 we obtain 



(55) 



Sit,Unit)) < 5(T", C/„(r")) + o(e-W"-"?-**). 



(56) 



Now notice that J2f=i [^iRjit)) + ^QiRjit)) + Vj ■ P(i?j(t))j is a time inde- 
pendent quantity. Therefore, 

N , 

J2 mm + -QiRjit)) + V, ■ PiRAt)) 

3 = 1 ^ 

N , . 
= E m^Tn)) + ^QiRATn)) + V, ■ PiR,iTn)) 

= 5(r„,c/„(r„)). 
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Combined with Lemma and (^6|) , this implies 
H(T„(i), T„(t)) = S{t, Unit)) - 5(T„, [/„(T„)) + o 



By Lemma |l^ we get 



Hence (pSf) is proved. Now we have 



\Un - R\\h1xL'^ < 2||^- i?|l^ixL2 + 2||T„||^ix^2, 



such that, by (|55|) and (B6|) we infer 



N 



Choosing t large enough we have 

This concludes the proof. □ 
6.2 Coercivity 

From now on and until the end of this paper, the subscript n is removed when 
there is no possible confusion. For example, J7„ is now denoted simply by U . 
We first prove Lemma O. 



Proof of Lemma li . From Lemma ^. we already know that for any j — 1, . . . , N 
we have 

(s"{R,{t))Tit),Tit))>C\\T{t)\\l,^^.. 

where the dependency of S in j is understood (recall that S = E + + Vj ■ P). 
We remark that 

'Q'^{Rj)T{t), T{t)) = Im f Ti{t)T2{t)^odx 

= Im / (V^Ti(t)(V^T2(t))da: = (Q"{R,)^jT{t), ^,T{t)) . 

.had \ ' 
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Similar computations can be performed for the momentum and the 0-order part 
of the energy. We deal with the gradient part by means of the classical IMS 
localization formula (see e.g. 



N , 
.7 = 1 ^ 



Straightforward computations using the definition of the cutoff functions <\)j 
and imply that 



< 



This implies that 



Combining these informations, we infer that 



Vt 



c 



N 



H{T{t), Tit)) = J2 {s'jmt))T{t),T{t) 



N 



J2 {s"iR,{t))^T{t), ^Tit)) - 



> c- 



Ti 



where the last inequality follows from the fact that t > Tq and Tq can be chosen 

2 



so that T( 



o>(^)' 



□ 



6.3 Almost conservation 

In this section, we prove Lemma R ecall that we have assumed that U = Un 
verifies the bootstrap assumption (|50|). We start with a preliminary lemma. 

Lemma 15. Let (j) he a function of the variable xi such that (f> and (f)' are 
bounded . Then for all t in the time interval of existence of U we have 



d_ 

dt 



Im / uiU24>{xi)dx = Its\ I dx^uiuicf)' {xi)dx^ 



— Re I dxjUiU24'ixi)dx ^ -Re I dx^uidx^uicj) {xi)dx 7^ 1), 



d_ 
dt 
d_ 
dt 



Re / dxiUiU24>{xi)dx 



p+l 
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Proof. The results follows from elementary computations using the fact that U 
is a solution to (0). □ 



Proof of Lemma 15. Let U — (ui, U2) and let us start by looking at the deriva- 
tive of the localized charge. By Lemma n5^ we have 



— Im / uiU2'4'idx 
at ' 



Im / dx^uiuidx^il'jdx / uiU2dtipjdx 

I f — rrijt 



< 



Im 


f 


1 






[m 




c 


L 


VtJ 


A, 



[dx^uiui —U1U2 ] V' 

(|9j;,Mlf2i| + \uiU2\)dx, 



dx 



where Aj := {x € R ; ^j{x) 7^ 0}. Remembering that (f>j = ipj — ipj+i for 
j = 1, . . . , iV — 1 and (f>N = ipN, and defining Aj := {x € M.'^; ip'j 7^ 0} we have 



dt 



Im / uiU2(f>jdx 



< 



71 



ildx^uiuil + \uiU2\)dx. 



and then 



dt 



Im / uiU24>jdx 



C ( 

< ^(ll«ill?/i(A,) + ll«2iii2(^^) 



Now notice that \\Utm(^A,)y.L^{A,) < '^\\U-R\\m(M'')y:mM^)+'^\\R\\m{Ai)xmA,)- 
Thanks to Lemma IS we have 



l^ll//i(Aj)xL2(Aj) - \\R3\\m(Aj)xL^{Aj) 



By using the properties of our partition of unity and the decay of the profile of 
Rj it follows that 



\Rj\\H^{Aj)xL'^{Aj) - 



Recall that a stems from (|52D. We conclude thanks to the bootstrap assump- 
tion @ that 

^Q,(t,C/W) ^ (57) 



if we choose a < Now we focus on the derivative of the localized momenta. 
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We start with the first component of the momentum. By Lemma |T5| we have 
■Re J dxiUiU2ipjdx 



d_ 

dt 



< 



Vt 



mi,^ 2 l"2| iVwip m|ui| 



|p+i 



< 



_c r 



P + i 

c r 



< 



Vt 



Now we argue as for the derivative of the locaHzed charge. The other components 
of the momentum can be estimated in a similar fashion. Thus we have 



Pj(t,U{t)) 



(58) 



Remark now that 



N 



S{t, U(t)) = EiUit)) + —Qj(.t, U{t)) + V, ■ P,{t, U{t)). 



Since is a conserved quantity, combining ( pTj ) and (^8|) gives the desired 
result. □ 



6.4 The Taylor expansion 

We now prove Lemma ^ We start by an estimate on the modulation parame- 
ters. 

Lemma 16. For any t e [i*,r"], we have 

N 



3 = i 



Proof. Recall that U = J2j=i -^i + Thanks to the interaction estimates given 
by Lemma |l^ and the orthogonality conditions it follows 

Q^{t, U) = Q{kj) + Im / TiT2^jdx + (^(e-^^V™— ^**). (59) 

We already computed the time-derivative of Qj during the proof of Lemma ^ 
(see (p7h), and it implies 



\Q,{t, U{t)) - Q,{r\ U{T^))\ = o(e-2-V^?^-**). 



(60) 
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Thanks to the scaling property (Q) of the profile we get 
Q{R,{t)) - Q{R,{T,,)) 

where the last inequality is due to the fact that u}j{T"-) — ojj. By simple Taylor 
expansion in frequencies we conclude 



2 d 



1] 



)(m-a;^)p-i 2- 



+ o(^i,(t)-c^,)- (61) 



Since cjj is part of the set Cstab (see ([l|)) we have 



> 0. 



Combining the bootstrap assumption (pd), and (p9|)-(pl|) gives the desired result 



□ 



Proof of Lemma 14- The first step consists in splitting the action using U 



+ ^- start with the energy part. We have 



N 



Y,E,{t,U)^E{U)^E{R + r} 



E(R) + E'{R)T + - (^E"{R)T, T) + o(|lT|l^,,^.)- 



We treat the order term first. By Lemma |T^, we have 
E{R) 



N 



2 21, 



2 p + 
1 



p+1 



p+1 
p+1 



where w^, = max{|(I)j|; j = 1, . . . , N}. In short, we have 

N 



37 



Now notice that 



(a;^-tJ^)+o(|c;}^-tj^|), 



such that, thanks to Lemma we get 

N 



E{R) = J2E{Rj) + Oie-^'^V"'-'^*''*'). 



Using similar arguments we have 

N 



i;"(i?)T,Tj) = ^(|£;"(i?j)T,T^ + 0(e-3"\A^^'**). 

The proof follows the same steps for the localized charges and momenta: we 
have 



N , . 



The second step consists in expanding ujj around ujj using Lemma 16 Re- 



membering that Rj is a critical point of E + ^-Q + ■ P, we infer 
E'{R,) + ^Q'iR,) + V, ■ P'{R,) = '^^1^^Q\R^). 



I3 



From Lemma 16; (|5C|) and (^4|), it follows that 

w-j- — Co 



7j 



Q\R,)T 



< 



The only thing left to see is to remove the tildes corresponding to modulation. 
We have 



N 



Y^{e{R,) + ^Q{R,) + v,-P,{R, 



N 



3 = i 
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where we have used the fact that 

{Cj, - u,){E' + ^g' + V, ■ P'){Rs)^ = 0. 

Thanks to Lemma |l^ we have 

N 

Gathering all these informations we get the desired result. □ 
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A Appendix 

Lemma 17 (Rellich-Kondrachov in H^). Let U be a bounded open set, s > 
and e > 0, and Un G H'^{W^) be a bounded sequence such that supp Un C il. 
Then there exists u G such that ||u„ — u\\^s-e = o(l) 

Proof. Let m„ be a bounded sequence in H^{D,) weakly converging to u G H", 
we shall prove that, up to subsequences. ||w„ — u\\jjs^^ = o(l). By Plancherel 
identity we have 



\i\<R 



+ {i + \ey-'\unio-m\'dc 

J\i\>R 



We have 

{i + \er-'\un{o-m\'d^< 



i?i>fl 

1 



(i + \er\un{o~um'd^< 



(i + R^y J ' ^ ' ^ - + 

in addition f2 is bounded and by weak convergence we have UniO ~^ ^iO- To 
conclude it suffices to show that 

{i + \er-'\i^n{o~u{o\'d^ = o{i). (62) 

i?i<fl. 

Notice that 

||Wn|lLoo(o) < l|Wn|lil(o) < K^)^ \\Un\\ l2 (^n) < K^)^\Wn\\H^ 
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and hence (1 + - "(OP is dominated by C(l + \R\'^y-^ such that 

(H) holds. 

Now, fix (5 > and choose R > and N sufliciently large such that 



iei>fl. 



and for all n > iV 



\i\<n 



l.e Un - u 



< d. 



□ 



Lemma 18 (Interactions estimates). There exists f e L^L;^(I 
such that if i ^ k 



\R,Rk\ + \Rj^Rk\ + |Vi?jVi?fe| + + |Vi?j|0fc < Ce'^''V^^^''^'f(t,x) 



N 



\R\P+^ -Y^\Rl\P+^ 



1=1 



N 



\R\p-'R-j2\Rir'Ri 



1=1 



N 



1=1 



Proo/. We start proving that there exists / G Lf LKR^R"^) n L^L^f (M,K'^) 
such that li i ^ k 

\RjRk\ < Ce"W"--'^**/(t,x). 

Thanks to (^ (the Lorenz transform gives indeed only a contraction along the 
direction of propagation) we know that 



\Rj\ < Ce" ^"""^3^'^"^^*' < C'e"5V™^"?l^-''j*l 
By a simple change of variable we get 



\Rj\\Rk\ < Ce-5V"-"2l=^le-3V"-"?l^-K-''')*)l, 
such that, thanks to the following inequality 

\X - {Vk ~ Vj)t\ > \{vk ~ Vj)t\ - \x\ > Vi,t - \x\ 

we conclude 

\Rj\\Rk\ < Ce^3V'»-'^?l^le-3V™-'^>**. 

Taking 3a < j we get the desired estimate. The estimates for |i?jVi?fc| and 
|Vi?jVi?fe| follow analogously. 
Now we shall prove that ii j ^ k 
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Let us suppose without any lack of generality that j < k — 1. Notice that 
that implies 

By a simple change of variable we get 
Now, for t > max{^, 1}, it follows 

< C'e~5V^™~^''**e~^V^™~^''^'. (63) 

Now for a < ^ we conclude 

The case j > k—1, j ^ k, follows identically as well as the estimates concerning 
the gradient. The second part of the lemma follows from the inequality 

i\a + b\P - \a\P - \b\P) < C{\a\\b\P-^ + |a|P->|) with p > 

that derives from the elementary inequality 

(|1 + t\P-l- \t\P) < C{\t\ + \t\P-^) with p>0. 

By arguing as before we get the desired estimates. □ 
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